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ABSTRACT 

The spectral and timing properties of an oscillating hot thermal corona are investi- 
gated. This oscillation is assumed to be due to a magneto-acoustic wave propagating 
within the corona and triggered by an external, non specified, excitation. A cylindrical 
geometry is adopted and, neglecting the rotation, the wave equation is solved in for 
different boundary conditions. The resulting X-ray luminosity, through thermal comp- 
tonization of embedded soft photons, is then computed, first analytically, assuming 
linear dependence between the local pressure disturbance and the radiative modu- 
lation. These calculations are also compared to Monte-Carlo simulations. The main 
results of this study are: (1) the corona plays the role of a low band-pass medium, 
its response to a white noise excitation being a flat top noise Power Spectral Den- 
sity (PSD) at low frequencies and a red noise at high frequency, (2) resonant peaks 
are present in the PSD. Their powers depend on the boundary conditions chosen 
and, more specifically, on the impedance adaptation with the external medium at the 
corona inner boundary. (3) The flat top noise level and break as well as the resonant 
peak frequencies are inversely proportional to the external radius r y (4) Computed 
rms and f-spectra exibit an overall increase of the variability with energy. Comparison 
with observed variability features, especially in the hard intermediate states of X-ray 
binaries are discussed. 
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1 INTRODUCTION 

X-ray binaries (XRBs) exhibit large variability on various 
timescales. While their spectral states and their accretion 
rates are typically changing from weeks to days (see e.g. the 
different c anonical states observed in black hole binaries, as 
defined bvlMcCIintock fc Remillardll2003l , iHoman fc Bellonil 
l2005l . lBellonill2009l and references therein), their light curves 



exhibit drastic changes from hours to millisecond. Several 
tools are now used in order to analyse these timing features 
(see e.g. Ivan der Klisl 1200*41 for a review) such as e.g. cross- 
correlation between different energy band (which allows to 
infer the so-called time-lags), but the most commonly used 
is still the Fourier analysis via the computation of Power 
Spectral Distribution (hereafter PSD). 

Despite its known limitations (e.g. signal phase lost in 
the analysis), any attempt to model the physics of XRBs 
has to take into account the various features observed in the 
PSD and its evolution when the source transits from a state 
to another. In the so-called hard state, the level of the vari- 
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ability is high and the PSD harbour a Band Limited Noises 
(BLNs) shape extending up to a break frequency Vb- Large 
peaks, the so called Quasi- Periodic Oscillations (QPO), can 
also be observed. In the soft state the Poissonian noise is 
usually dominating on all the frequency range and no or 
weak QPOs are detected. Different types of QPOs (called 
A, B and C) can be identified depending on the value of 
their frequency, their st re ngth and even their time lags (see 
again Ivan der Klisll2004l or lCasella et al. I l2004l and references 
therein). However, if QPOs are remarkable features, the 
major part of the variability is usually aperiodic. As for en- 
ergy spectra, variability evolves during time: in general, the 
overall variability decreases as the spectrum softens, with 
frequencies increasing, unt il in the soft state where the BLN 
reaches a low level (see e.g. lBelloni et aI.ll2005l ; lBeIlonill2009h . 



Several models intending to interpret variability 
features focuses on the QPO phenomena. For the high 
frequen cy QPOs (v > 100 Hz), lense-thirring precession 
(see e.g. lStella fc Vietril 19981) or beati ng frequency between 
particular orbits ( Lamb fc Millerll2003l l have been proposed. 
For the low frequency QPOs (0.01 < v < 100 Hz), fewer 
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Figure 1. Sketch of the model geometry. The corona is assumed to harbour cylindrical symmetry. It extends from its inner radius ri to 
its outer radius r; and its height is h c . Non-specified white noise excitations are assumed to occur at r; then triggering (magneto-)sonic 
waves in the corona. The resulting X-ray luminosity is then computed through thermal comptonization of embedded soft photons. 



models are available (see e.g. ITagger fc PellatI 1 19991 or 
Titarchuk fc Shaposhnikovl 120051 . but also IStella fc Vietril 



19981 ). but they usually try to explain the observed correla- 



tion between frequencies and/or other observables, without 
taking into account the whole emitting process. However, 
it is worth noting that it is in the highest energy bands 
that the X-ray flux is observed as being highly variable. In 
contrast in most of the available models, the source of vari- 
ability lies in the geometrically thin accretion disc, which 
emit mainly at lower energies. Note that some observations 
also suggest dir ectly that the disc is less variable than the 
corona (see e.g. IChurazov et al.ll200ll , lRodriguez et alj|2003l 
however see also Wilkinson fc Uttlevll2009i r In consequence 



a proper model for variability in BHB has to deal with the 
radiative transfer between the disc and the corona. 



A few models have already been proposed in the liter- 
ature. For instance, a fu ll mode l ling o f the BLN component 
has been attempted by [Misra ( 2 000l) in the fram e work of 
a "transition disc model": following iNowak et ail |l999h 's 
idea, an acoustic wave is propagating within the accretion 
disc. In this framework, the author considers only the propa- 
gation to occur in one direction towards a central sink. This 
model is then used to explain qualitatively the general shape 
of the power spectra and lag energy dependencies observed 
in Cyg X-1. 

IPsaltis fc Normanl |2000l ) tried to model the filter effect of 
a narrow annulus in a geometrically thin disc. However, the 
nature of this annulus is not specified and its typical exten- 
sion is Sr/r < 10 -2 . They also neglect the contribution of 
the radial pressure forces. In this framework and depending 
on the mode of oscillation chosen, the external perturbation 
is shown to be modulated in amplitude according to the 
exciting frequency. The square of the response in pressure 
exhibits Lorentzians which could account for the observed 
behav iour in BH or NS binaries. In a similar spirit. lLee et al.l 
|2004l ') studied numerically the response to radial perturba- 
tions in an accretion torus. They showed that resonances 
could occur and become larger when the frequency differ- 



ence between the radial and vertical epicyclic frequencies 
where half the forcing frequencies. 

One of the very first attempt to evaluate the effect of 
propagation in a comptoni sing region on t h e tim ing be- 
haviour was examined by IMivamoto et al.l (|l988l ). They 
showed that the tight period dependence of the observed 
time lags in Cyg X-1 could not be accounted for by the 
inver se Compton scatterings proc ess only. In another frame- 
work, [Zvckr&_SobcfcwsJ{g (|2005h tested their model of vari- 
ability (jzvckill2003f r where the variable emission responsi- 
ble for the noise component is attributed to multiple active 
regions/perturbations moving radially towards the central 
black hole. The QPOs are obtained by modulating either 
the reflection amplitude, the heating rate, the covering fac- 
tor of the reprocessor or the column density value. They 
predict for each of these cases the corresp onding power and 
f spe ctra (for a definition of f-spectra, see IRevnivtsev et akl 
Il999l ). and t he time lag energy depend encies. In a subse- 
quent paper. ISobolewska fc Zvckil |2006l ) tested their model 
on real data and concluded that within this framework the 
QPO spectra in the hard state are always softer than nor- 
mal average spectra and cold disc oscillations might then be 
responsible for the low frequency QPO. On the other hand, 
when the energey spectra gets softer, the QPO spectra are 
harder and the low frequency QPO might then originate 
from the hot plasma. 

In a more recent study, ISchnittman et al] l|2006l ) try to 
model the oscillation of a torus in Kerr metric, via three 
radiative processes: a thin emission line then a thick one, 
and finally an optically thick thermal emission process. They 
manage to reproduce some of the the properties of C-type 
QPOs, especially the observed increase of the amplitude 
with the inclination of the system. 

In this paper we present a new approach that deals with 
the radiative transfer into an oscillating corona. These oscil- 
lations are assumed to be due to a magneto-acoustic wave 
propagating within the corona, modulating the efficiency 
of the comptonisation process on embedded soft photons. 
These basic ingredients give a promising framework to re- 
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produce the main timing features of the X-ray binary in 
hard and hard-intermediate states such as the Band Limited 
Noise continuum and C-type LFQPO. The assumptions of 
the models and wave equation solutions are detailed in sec- 
tion [2] An analytical study in the linear approximation is 
discussed in section[3]and compared with Monte-Carlo sim- 
ulations in section [4] We finally discuss the main results of 
this toy model and the comparison to observation in section 

El 



2 THE MODEL 



2.1 Basic assumptions 



The general structure of the model is sketched in Fig. [T] A 
hot optically thin cylindrical medium is assumed, hereafter 
called the "corona". It is limited radially by its inner and 
outer radii n and rj respectively. This corona has thus a 
ring shape of typical height h c . This geometry agrees with 
the observations that suggest the X-ray corona to be closely 
linked to the base of the jet in XRB (e.g . iMarkoff et~alll200ll ; 
iFender et al. 1 12004 IMarkoff et al.ll2005T ). Then the inner ra- 
dius n can be identified as the last stable orbit close to the 
central compact object. On the other hand rj can be com- 
pared to the transition radius between an outer standard 
accretion disc and the inner hot corona, a geometry com- 
monly invoked in XRB (e.g. lEsin et ai]|l997l . I Ferreira et al.l 
120061 . iDone et al.ll2007h . 

In this paper, no hypothesis is done on the physical ori- 
gin of these radii. For simplicity, the corona is assumed to 
have, at rest, a constant temperature To and density no and 
consequently a constant pressure Po- We then consider that 
pressure instabilities at the external radius rj generate a 
sound wave within this thermalised plasma. The origin and 
the nature of those instabilities are not discussed in this 
paper since we focus only on the radiative response of the 
corona. We assume these instabilities to have a white noise 
spectrum (i.e. same amplitude for all excitation frequencies) 
which corresponds to a Dirac perturbation in the temporal 
domain, i.e, we limit our study to the corona transfert func- 
tion . For seek of simplicity we also restrict our calculation 
to the ID case i.e. the wave will only propagate radially at 
the sound velocity: 




3.1 x 10 s 



To 



100 keV 



1/2 



(1) 



In the last expression To is in keV. It is important to note 
that there are simplifications which make this paper only a 
first step. Once additional complications such as incorporat- 
ing the effects of rotation are introduced, they might lead 
to some revision of the scenario presented here. This will be 
addressed in a forthcoming paper. Note however that the 
nature of the corona and hence its exact rotation profile is 
still mainly unknown. 

Finally we suppose a blackbody seed photon field of temper- 
ature Tseed emitted isotropically at the corona midplane (i.e. 
in z = see Fig.[T]). These seed photons will be comptonized 
in the corona then producing a variable X-ray emission. 



2.2 Wave equation 

In cylindrical ID geometry, in absence of local damping or 
excitation, and neglecting the rotation, the basic wave prop- 
agation equation can be written as: 



r dr I dr J c? dt 2 



(2) 



p(r, i) is the perturbation in pressure given by, in complex 
notation and for a given frequency v — uj/2n: 



p(r, t) — p r (r) exp 1 ' wt ' . 



(3) 



Introducing the new variable x = — r and putting Eq. [3] in 

i-i Cs 

Eq. [2] we get: 



,d p r dp r 2 

-3—5- + x— t- x p r 

ax 2 ax 







(4) 



The general solution of those equations are linear combina- 
tion of Hankel's function of the zeroth order Hq and Hq 
jAbramowitz fc Stegunlll964l 'l: 



p(x,t)= [aHl{x) + f3Hl{x)} e 



(5) 



a and /3 being determined by the boundary conditions. Note 
that the above sol utio n is similar to the one obtain by 
iNowak et all ij 19991 ). or iTitarchuk fe Shaposhnikovl ()2005l ) 
despite the fact that those authors used only one instead of 
a linear combination of both Hankel's function as solution. 



2.3 Boundary conditions: total reflection at the 
internal radius 

At the external radius of the corona rj , we assume a constant 
excitation pJ3 (white noise hypothesis): 



— ZU)t —lUlt 

e = poe 



(6) 



where Xj = x(r = 



The closure relationship will be 



given by the behaviour of the wave at the internal radius 
n. We will, as a primary assumption, consider that there is 
no transmission of the wave and hence total reflection in n. 
If n is equal or close to the Last Stable Orbit, one would ex- 
pect indeed the density of the corona to drop quickly inside 
this radius and hence the pressure as well. It thus gives: 



p(xi,t) = (aHl(x % ) + /3Hi(x t )) e 



0. 



where Xi — x(r = r;).The system of Eqs. [S] and [7] can be 
solved in order to obtain a and (3: 

Hl{xi) 



Hlix^HUxi) - Hl{ Xi )H*{ Xj ) 



(8) 
(0) 



3 ANALYTICAL SOLUTIONS IN THE LINEAR 
AND ZERO-PHASE APPROXIMATIONS 

In this section, we infer the shape of the expected PSD in 
the simple case where the radiative response of the corona 
depends linearly on the local perturbation. This appears to 

Note the difference between po, the pressure perturbation im- 
posed in r = rj and Pq the corona pressure at rest 



4 C. Cabanac, G. Henri, P.-O. Petrucci, J. Malzac, J. Ferreira, T. M. Belloni 




0.001 0.01 0.1 1.0 10.0 100.0 1000. 0.001 0.01 0.1 1.0 10.0 100.0 1000. 0.001 0.01 0.1 1.0 10.0 100.0 1000. 



v/Vj V/Vj v (Hz) 

Figure 2. Left: Power spectrum plotted in reduced frequency Xj = v/vj (with Uj = 2nc 3 /rj) for a single zone corona with £ = 10 — 2 , 
Clin = 1 an d SPq/Pq = 0.1. Center: same spectrum plotted in vP v . Right: PSD plotted in true frequency and its evolution when 
moving rj inwards only. Here, M = 10 M Q , n = 2 R g = 3 X 10 6 cm and kTo = 100 keV (hence c s = 3.1 X 10 8 cm.s -1 ). Black: 
r i = 680 R g = 10 9 cm, red: rj = 200 R g = 3 X 10 8 cm, blue: r-j = 68 R g = 10 8 cm, green: rj = 14 R g = 2 X 10 7 cm. (Note that 
(%rms) 2 (l v = 0) slightly decreases and fb increases when the external radius of the corona decreases, though the effect is tiny here.). 



give analytical results in good agreement with our Monte 
Carlo simulations detailed in the next section. 

The corona being optically thin, the luminosity dLo 
emitted locally at rest in a ring of radius r, width dr and 
height h c is proportional to its emissivity per unit volume 
Vo- 



with L = m-] h c (rf — rf) = Ttrfoh c x 2 j^{l - and 



dLo(r) = rjo2irrdrh c 



(10) 



In presence of the sonic wave, the luminosity dL(r,t) of this 
ring varies in time and we will make the assumption that its 
relative variation (dL — dLo)/dLo is a linear function of the 
relative pressure perturbation p(r,t)/Po: 



dL(r,t) - dL (r) p(r,t) 

— L-Hn ~ • 

-TO 



dLo(r) 

din is the constant of proportionality which is in complete 
generality a function of the perturbation frequency. How- 
ever, if the delays implied by the multiple diffusions of the 
photon inside the corona are low compared to the period of 
the wave, little phase delay is expected between the pres- 
sure wave and the luminosity. Hence Cu n will be a real 
number independent of oj. This is what we call the "zero- 
phase approximation" . It is usually verified in optically thin 
pla sma where photons travel a few h c before escaping (see 
e.g. iMalzac fc JourdainlfeoOOl) . Then a rough estimate of the 
time spent by the photons inside the corona compared to 
the period 2n/oj of the wave gives: 

uj h c _ 1 h c c s 

2tt c 2-7T r c 
It is generally much smaller than 1 in the cases we are 
interested in (i.e. corona aspect ratio h c /r < 1 and corona 
temperature of a few tens to hundred of keV) unless x 
becomes of the order of 100 or 1000. 



(12) 



Combining Eqs llOl and llll we obtain: 

dL(r,t) = V o (l + CuJ^- ] 



(13) 



Previous equation integrated on the whole volume of the 
corona gives: 

2-Kh c Clin r j 



L = L + 770 - 



Po 



p(r, t)rdr, 



(14) 



? = - = (15) 

Note that relationship [TJ is only valid if the emissivity per 
unit volume 770 of the corona, at rest, is uniform, which 
is a direct consequence of our assumptions of constant 
temperature and density Tb and no. 



3.1 Power Spectra obtained with large vertical 
wavelength and total reflection in n 

Using Eq.0 the expression of the perturbation in luminosity 
L~ becomes 



L, 



2LoCnn e 



Po 



(aHl (x) + (x)) xdx. 



(16) 

The Hankel's functions pr esent in a and /3 (see Eg s. [5] and O 
are easily integrated (see lAbramowitz fc Stegun|]l964l ') and 
hence the previous equation leads to: 



2LoCunPo e 



x 2 Ai~C< 



-M. 



with, 



M x s = x 



rr2 rrl rrl rr2 



(17) 



(18) 



H\ and H\ are the Hankel's function of first order and the 
index i or j corresponds to the point where the Hankel's 
function is evaluated i.e. Xi or Xj respectively. It is also 
easy to demonstrate that the function M is real in case of 
total reflection (both numerator and denominator are pure 
imaginaries and hence the ratio is real). 

The PSD P v = \L^\ 2 /L can then be directly deduced 
from Eq. [T7] 



P v = \L^\ 2 /Lq — \2Ciin 



po M(x s ,0 



P x 2 Jl-e 



(19) 



It is directly proportional, due to our linear assumption, to 
the input perturbation through the term 2Cu n po / Po but it 
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Figure 3. Dependency of the reduced PSD power at lowest fre- 
Pxj (xj ->• 0) 

quencies (xj — > 0) — — ^ (see Eq. 1 1 9 6 in function of the 



ratio £ = fijrj 



Clin p 



is also modulated by the intrinsic response of the corona 
through the function M(xj,£)/(x*(l - f; 2 )). 

In the zero-phase approximation (see above) Cu n is real 
and we can find asymptotic expressions of the PSD for low 
and high frequencies. For low frequencies, the function M 
can be approximated to (see appendix [A}: 



(20) 



Hence, the PSD P„ (which is inversely proportional to x 2 ) 
tends to a constant. At high frequencies, the Hankel's func- 
tions tend to cosine functions whose amplitudes are propor- 
tional to xj 1 ^ 2 and thus P v oc xj 2 . 

Examples of PSD given by Eq. \T§\ are plotted in Fig. 
[2] in function of the reduced frequency Xj — vjvj, with 
Vj — 2nc s /rj. As expected a flat-top noise component at 
low frequencies and a red noise at high frequencies are 
present. The break frequency scales like Ub ~ 2.5i^ oc c s /rj 
when £ — ¥ and v\, ~ nvj/(l — £) when £ — > 1 (see some 
examples plotted in right panel of Fig. [2]). To determine this 
value of the break frequency, we fitted the analytical PSD 
obtained with a zero c e ntere d Lorentzian, following the defi- 
nition of lBelloni et al.l (|2002l '). The PSD also exhibits a peak 
around the break frequency and oscillations above it as ex- 
pected from Hankel's function. The ratio between the fre- 
quencies of the first peak and the second one is predicted 
to be close to 2. Note therefore that the frequencies depend 
on the value of the outer radius of the inner region. We also 
noticed that in our model the QPO peaks around the break 
frequency. 

Moreover, the PSD plateau, at low frequencies, slightly 
decreases when £ increases. This is due to the behavior of the 
modulation function M(xj, £,)/(x 2 (l — £ 2 )) at low frequency 
(i.e. when Xj tends to 0). It is plotted versus £ for Xj = 
in Fig. [3] and shows a decrease from about 25% for £ = to 
about 5% for £ = 1. 



3.2 Wave transmission in r\ 

The oscillations present above the break frequency in the 
PSD (see Fig. are due to the infinite resonances at the 
cigenfrequencies of the corona due to the assumption of total 
reflection in n and no wave damping. This total reflection 
hypothesis can be relaxed by assuming that part of the wave 
is transmitted in n in a medium of different sound velocity. 
An impedance adaptation then occurs between the two me- 
dia, depending on the value of the acoustic impedance of 
the system Z = c 3 ,2/c 3 ,i, i.e. the ratio of the corona sound 
speed c s ,i to the sound speed c s , 2 below n . The general 
form of the solution in the corona (medium 1) is similar to 
the one obtained before (Eq. [5]) and hence can be written in 
reduced units x = ujr/c Si \ as: 

pi(x,t) = [aH% (x) + ml (*)] e~ mt , (21) 

whereas, in the medium 2, the transmitted wave is only pro- 
gressive and hence has the following form: 



Pi(x,t) = (Ho e" 



(22) 



In order to constrain a, /3 and (, we need now three different 
equations. The first ones comes from the pressure continuity 
in Ti and r,: 



aH^Xi) + /3H%(xi) = C#o(J),and (23) 

aH^x^ + pHoixj) = Po- (24) 
The mass conservations in r; gives the third relation: 

(/OlSlUl),., = (p2S2V 2 ) ri , (25) 

with 5*i = 2irrih\ (respectively S2 = l-nrih^) being the 
vertical surface in n in medium 1 (respectively medium 
2) and Vi (respectively V2) the corresponding flow velocity. 
The link between v m and the pressure perturbations p m in 
each medium is obtained by applying the Euler equations 
(m = {l,2}) : 

P m ~KT = — VjJm- (26) 



dt 



Hence, by using Eqs. I21l and 1221 we get: 

[aHl(x) + f3H 2 (x)] and 

7 [<-*©]■ 



ujpiVi 



-iujp 2 v 2 — 



Cs.l 
UJ 

ZCs 



Combining these two equations with Eq. [25] give then: 

"Si r rrlr \ , ott2, \1 ^2 



\aH{{xi) + pHl(xi 



(27) 
(28) 

(29) 



C s ,l L " Zc s 

We need then to constrain the corona height in both side 
of ri. We therefore assume that the corona is in hydrostatic 
equilibrium and in Keplerian motion. Consequently, h/r = 
Cs/vKepiir) in each medium. As a result, at the internal 
radius n, S m /c s , m = 2nnh m /c s ,m = 2tyt 2 /v K e P i(ri). This 
latter value is independent from the value of m and hence 
we obtain that Si/c B: i = 5*2/c s> 2 = S2/(Zc B: i). The mass 
conservation (Eq. I29|l therefore reduces to: 

aHl ( Xi ) + /3H?(x t ) = CB\ (|) , (30) 

The resolution of the previous system gives then the full 
solution in pressure within the corona: 

AHl{x) + BH%{x) 



Pi{x,t) = po 



(31) 
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Figure 4. Left: power spectrum obtained for the corona when the sound velocities in media 1 and 2 are the same (Z = c Sj 2/c Sj i = 1). 
Right: PSD evolution when changing the value of the acoustic impedance. Thick grey curve: Z = 1.2, black: Z = 10. 



with, 



A = Hl(t-Xi)Hl ( fy) -HlUx^Hl ( fy) and, (32) 



B = Hl^jHfaxA-Ht^jHlitxj). (33) 

In the linear approximation, we can obtain the outing PSD 
by following the same steps as in section 13.11 The power 
spectrum from the main corona is therefore very similar to 
Eq.[l9] 



p x 2 Ji-e 



with 

C 
D 



C + D 
Xj AH&xs) + BH*( Xj y 
A[Hl( Xj ) - tHlfaj)] and 
B [H!( Xj ) - ZHlitxj)] , 



(34) 



(35) 

(36) 
(37) 



A and B being given by Eqs. [32] and 1331 

Some examples of PSD are plotted in Fig.Ufor different 
values of the acoustic impedance Z — c s ,2/c s ,i. As expected 
the resonances are cancelled when there is impedance adap- 
tation (i.e., Z — 1) between both media (see left panel of 
Fig. 2}. The sound wave is in that case totally transmitted 
in n. Then the strength of the resonances are tuned by the 
impedance value as soon as it is different from unity and 
the higher the difference between both sound velocity, the 
stronger the resonances. Two PSD examples with Z = 1.2 
and Z = 10 are displayed on the right panel of Fig. \$$. 



the corona mid plane [z = 0). The temperature and optical 
depth of the corona at rest are fixed to To and to respectively 
implying a uniform pressure Po- The wave propagation be- 
ing supposed adiabatic, the perturbation in optical depth St 
and temperature ST are given by: 

ST 7-1 p 

% ~ 7 Po' 
St 1 p 

to 7 Po ' 



and 



(38) 
(39) 



Then we impose sine perturbations in r — Vj with dif- 
ferent frequencies but same relative amplitudes e = po/Po 
(white noise). We assume total reflection at the internal ra- 
dius and the pressure (and hence temperature and density) 
profile that the photon encounters while travelling within 
the corona is then given by Eqs. [5] [8] and [9] For those simu- 
lations, both the linear and zero-phase approximation stud- 
ied in the previous section were released. Hence the pho- 
ton will encounter during its travel, and scatterings after 
scatterings, either some part of the corona where the per- 
turbation in pressure is positive, or negative. The number 
of positive and negative zones in the corona is linked to 
the wavelength and hence the exciting frequency. The full 
Klein-Nishina cross section is taken into account. The tem- 
poral evolution of the corona during the photon motion is 
also fully taken into account but appears to be negligible 
compared to the wave frequencies used, in agreement with 
our zero-phase approximation adopted in Sect. [3] We finally 
fit the corona X-ray emission with a sine function in order 
to obtain its amplitude and potential phase delay. We 
repeat this procedure for different frequencies in order to 
build a power spectrum. 



4 MONTE-CARLO SIMULATIONS 

4.1 A linear comptonisation code 

In order to check the validity domain of the linear approx- 
imation that is used in the previous section, we performed 
Monte-Carlo simulations of the radiative response of the 
corona. The comptonisation code used is linear, i.e, no feed- 
back of the computed high energy flux on the state of the 
corona is taken into account. The code uses the weigh ted 
Monte-Carlo technique (see e.g. iPozdniakov et aL I ll983h . A 
monothermal blackbody distribution of temperature T seec i is 
assumed for the seed photons which are randomly drawn in 



4.2 Results 

Some examples of power spectra are plotted in Fig. [5] 
for different values of the external radius rj and different 
energy band indicated on the different figures. The other 
parameters are n = 3 x 10 6 cm, kTo = 100 keV, to = 1.4 
and kT scc d = 0.25 keV. The amplitude of the modulation in 
pressure is set to an arbitrary value e = 0.15. 

For such low value of e, the PSD behaviour appears in 
good agreement with the linear hypothesis and zero-phase 
approximation studied in the previous section. Note that, 
thanks to the Monte-Carlo simulations, we are also able to 
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V (Hz) v (Hz) - - v (Hz) 

Figure 5. Evolution of the power spectrum obtained with either rj or the energy band. 1 cr errors arising from the fit on the amplitude 
of the wave in flux are also plotted. See text for the input parameters value. Upper panel, from left to right: PSD obtained for 25-45 keV 
band. The QPO frequency moves upward from 0.33 to 0.67 and 3.3 Hz when the external radius of the corona rj moves inward. Lower 
panel: Evolution of the PSD with the energy band. The value of the external radius is r 3 ■ = 3 X 10 8 cm. Note the increase of the overall 
variability with the energy. In Fig. \5\ f The dashed curves are Lorentzians usually employed to fit the different components in PSDs of 
XRB. 




Energy (keV) Energy (keV) Energy (keV) 



Figure 6. Left: Cu n (circle) vs C qua d (cross) value with linear coordinates on the y axis (see text for the parameter of the corona in the 
steady state). Under about 3 keV (the "pivot") the value of Cu„ is negativeThis is a direct consequence of photon number conservation 
in the compton scattering (see text). Center: same as right plot, but in absolute value and y axis in log coordinates, to compare the 
magnitude. Above 10 keV, Ca„ is higher than 1, emphasising the effect of modulation in the flux by the perturbation in pressure. Right: 
Plot of 

tmax spectrum (see text for explanation). Here, \th — 0.2. If e.g. e — 0.15 as plotted (dashed horizontal line), two linear and non 
linear zones can be observed. 
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study the evolution of the power spectra with energy. In- 
terestingly, the overall shape of the power spectra can also 
be mimicked by using sets of Lorentzians especially at high 
energy as shown in Fig. [5] f . 



4.3 When does the linear approximation become 
invalid? 

As shown in the previous section, for low values of the mod- 
ulation amplitude e, the Compton emission of the corona 
agrees relatively well with the linear approximation. We ex- 
pect however some deviation from linearity for larger mod- 
ulation amplitudes, deviation that should also depends on 
the energy. This aspect is investigated here by increasing e 
in the simulations and then adjusting the corona spectral 
emission at different energy E with a quadratic polynomial, 
i.e.: 



To further investigate the domains of non linearities, 
let's assume an arbitrary threshold \th of the ratio of the 



L(e,t)-L 



L 



8L(e,t) 

L E 
[Cun(E)e + C qua d{E)e 2 ] e 



(40) 



Hence, the energy-dependent parameters Cu n vs C qua d 
allows to "quantify" the linear and quadratic behaviour of 
the perturbation^]. 

Starting from a reference spectrum corresponding to 
fcTseed = 0.75 keV, kTo = 75 keV and r = 1.2, we simulate 
9 different spectra with values of e ranging from to 60% 
and fit the data with the polynomial given by Eq. [40] for 
different energy bins. The corresponding values of Cu n and 
C qua d are plotted in Fig. [6] in function of the energy. They 
appear strongly energy dependent. Noticeably they both 
cancel and change sign at a medium energy of about 2 keV, 
signature of a "pivot" in the variable spectral emission of 
the corona. 

This pivot is also clearly visible when comparing light 
curves at low and high energies such as those plotted in 
Fig. [7] Whereas the high energy lightcurve respond to the 
excitation coherently, the low energy one (under the pivot) 
has a phase lag of n. This pivot is due to the compton up- 
scattering of low energy photons which naturally produces 
a decrease of the number of soft photons and, simultane- 
ously, an increase of the number of the high energy ones. 
Consequently there is a n phase lag between the two energy 
domains . The quadratic term C qua d cancels a second time at 
higher energy (~ 90 keV in the simulation plotted in Fig. [6]). 

By examining the trend at high energy in center panel 
of Fig. [6l it is worth noting that the value of Cu n is higher 
than 1 above 10 keV. This means that any fluctuation in 
pressure will produce an amplified radiative response of 
the corona at high energy. This is due to the pivoting of 
the spectrum, which acts as a "lever" arm: the higher the 
energy bin, with respect to the pivot energy, the larger the 
luminosity variation. 



quadratic to the linear terms in Eq. 1401 i.e. 



,ad£ 



Clin 



above 



which we estimate that the linearity hypothesis is no more 
valid. This, in turn, determines different area in the e/xth 
vs energy plane (see Fig. [6] ) that agrees or disagrees with 
the linear approximation. For example, assuming a thresh- 
old \th of 20% means that we estimate that linearity is well 
satisfied as soon as ^- ade < g 2 . For the set of parameters 

Clin 

used in Sect. 14.21 and for which we obtain Fig. \6\ this con- 
strains the perturbation relative amplitude e to be smaller 
than £ mal ~ 50% = 2.5xth in the 0.1—3 keV energy band 
and even much smaller above 10 keV. For high values of e, 
say e.g. e = 0.7, the whole energy domain is non-linear. For 
lower values, say e.g. e = 0.15 as plotted by the horizontal 
line in right Fig. [6l two linear and two non linear zones can 
be noted. 



4.4 Energy dependence of the variability 

Finally, the strength of the variability and its evolution 
with energy is investigated. For that purpose, we com- 
puted and plotted on Fig. the so called "RMS spectra" 
and "f-spectra" . Those "variability spectra" were com- 
puted for three spectral ranges: at the lowest frequencies, 
where the PDS is flat (10~ 3 — 10 _1 Hz), near the resonance 
(10 _1 - 0.7 Hz) and after (0.7 - 1.5 Hz). 

As expected from the behaviour of Cu n , the presence 
of the pivot lead to a drastic decrease of the rms at around 
2 keV for all the frequencies studied, which is a strong pre- 
diction of the model, even if the position of this pivot in 
energy depends on several input parameter such as the seed 
photon temperature. After it, the variability increases with- 
out any cut-off at high energy. This due to the fact that the 
wave in the corona also modulates the temperature of the 
high energy electrons and hence the position of the cut-off 
in the high energy spectrum. The higher the modulation, 
the larger the cut-off. An increase of the variability of the 
0.7 — 1.5 Hz component is also predicted above ~ 100 keV. 

F-spectra were also computed, following definition of 
iRevnivtsev et af1 (|l999l ). and adopting the same three fre- 
quency bands as above. Such spectra were then fitted be- 
tween 2 and 150 keV by a powerlaw. Their photon indices 
are decreasing with increasing frequencies as they are equal 
to 1.61±0.02 for the 10" 3 - 10" 1 Hz band, 1.53±0.01 for the 
lO -1 -0.7 Hz band and 1.48±0.05 for the 0.7-1.5 Hz band. 
In contrast, the time averaged spectrum has a higher photon 
index (F a „ g ~ 2.1). The ratio of the obtained f-spectra to a 
powerlaw of index 1.53 is plotted in Fig. [7] and also the ratio 
of the time averaged spectrum to a powerlaw of index 2.1. 
The time averaged and the f-spectra exhibit a quite sensible 
deviation from the simple power-law model. This discrep- 
ancy with a power-law is however more emphasised in the 
f-spectra than in the time averaged spectra. 



2 Note that the parameter Cn n of Eq. 1401 is not exactly the 
same as the one used in Sect. [3] It also includes all the geometry- 

dependent part of the corona response (i.e. the term — = 1 ' 

present in Eq. 1171) 



5 DISCUSSION 

The timing response of an oscillating corona in cylindrical 
geometry is investigated, the corona acting as a filter when 
responding to a perturbation in pressure generated at 
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500 1000 1500 2000 0.1 1.0 10.0 100.0 0.1 1.0 10.0 100.0 



Time (s) Energy (keV) Energy (keV) 

Figure 7. Left: Light curves obtained for the lowest frequencies probed (y = 10 -3 Hz), for two energy bands, and their best fit by a 
sine function (dashed curve). Upper plot: 0.25 — 1 keV. Lower plot: 2 — 3 keV. Note the phase lag of 7r between the high and low energies. 
Center: "RMS spectra" for different frequency bands. ▲: 10~ 3 — lO^ 1 Hz, •: 10 — 1 — 0.7 Hz. □: 0.7— f.5 Hz. Right: 3 lower curves: Ratio 
of the f-spectra to a power law model with photon index a = 1.54, for different frequencies (same symbol code as center panel). Upper 
curve (*): ratio of the time averaged spectrum to a power law of index V a vg = 2.1. Time averaged and 0.7—1.5 Hz spectra were rescaled 
for clarity. 



its external radius. In the case of total reflection at the 
internal radius, infinite resonances may occur. They may 
be damped in the case of acoustic impedance adaptation 
between the corona and the external medium. We have 
shown that the filtering effect of the corona leads to power 
spectra with a broad band shape well fitted by a Lorentzian 
centered in v = (see Fig. [5}. Moreover, the resonances 
produce a major peak at the PSD break and smaller peaks 
at higher frequency. The characteristic frequencies of the 
PSD break and peaks scale with Vj = 2nc s /rj. For a 10 
solar mass black hole, an outer corona radius r,- = 200 R g 
and a corona temperature fcTrj = 100 keV (i.e. a sonic 
velocity c s = 3.1 x 10 s cm.s -1 ), this gives Vj = 0.17 Hz. 
These power spectra are thus relatively similar to the one 
observed in the HIMS (Hard Intermediate State) of X-ray 
transients, whose shape are well fitted by 3-4 Lorentzians 
peaking at different frequencies from a tenth of hertz (for 
the break and the low frequency QPO, hereafter LFQPO) 
to hundreds of hertz (for the high frequ ency QPO, see e.g. 
lNowak|[2000l or iPottschmidt et al J 12003 ). 
Interestingly our low frequency peak is always of the order 
of the frequency break Vb (both scaling with Vj ) . This is also 
in agreement with the observed co rrelations between ft and 
the l o w frequency QPO vlfqpo dWiinands fc van der Klisl 
1 19991 ; iBelloni et af]|2002l ; iKlein-Wolt fc van der Klidl2008h . 
We note however th at vlfqpo is usually 5 t i mes h igher 
than i/j,, especially in lWiinands fc van der Klisl (|1999| ) rela- 
tionship. We note ho wever that in the equ ivalent correlation 
plotted on Fig. 11 of lBelloni et al.l |2002l ). the lower branch 
shows a correlation between the break frequency v\, and 
the typical frequency of the "hump" (Ph), with a ratio 
close to one. Finally, the ratio between the first two peak 
frequencies is close to 2 in our model, as it is observed for 
several type C QPO. 



Our model also predicts a decreasing PSD power and 
an increasing resonance peak frequencies when decreasing 
the outer corona radius rj. This is here again consistent 
with the general trend observed in BHB where the PSD 
power decreases when the observed frequen cies increase (BH 
anti-correlation, IBelloni fc Hasingen 1 ll990T ) . However in our 
model, this effect is low and could hence only account for 



a fraction of the observed anti-correlation. Indeed, the BH 
anti-correlation translates the fact that in vP v , the different 
observed Lorentzians are peaking at roughly the same level, 
or equivalently that the rms integrated variability remains 
roughly constant. An extra condition on the exciting process 
is therefore required, with a rms level that must decrease 
when the size of the corona decrease as well. The nec- 
essary condition would hence be that SP oc rf , with a ^ 0.5. 

As shown in Fig. we expect a n phase lag between 
the very low energy and the high energy light curves, which 
is a strong prediction of the model. This is a direct conse- 
quence of the assumed geometry, especially of the fact that 
the source of seed soft photons is assumed at the corona 
midplane. This source is necessarily on the line of sight of 
the observer. Then due to the photon number conservation 
during the Compton scattering process, the disappearance 
of the soft photons is directly compensated by an increase 
of the hard ones. Moreover it depends on the pivot energy 
which crucially depends on the input parameters (especially 
feT c and kT so ft). In a geometry where the soft seed photons 
are produced outside the corona, this constraint can be eas- 
ily relaxed. This would be the case of a soft photon field 
produced by an outer accretion disc. Such geometry (inner 
hot corona surrounding by an outer accretion disc) is in- 
deed generally believed to qualitatively well reproduced the 
inn er region of the accretion flow aro und compact objects 
(i.e. lEsin et"alll 19971 ; Ibone et al.ll2007l ). Note that with such 
geometry, a decrease of the outer corona radius rj (which 
corresponds to the inner accretion disc one) would imply an 
increase of the seed photons temperature and flux i.e. an 
increase of the corona coolings. Consequently a softening of 
the X-ray spectrum should happen. Since a decrease of rj 
implies also, in our model, an increase of the resonance peak 
frequencies, a correlation between the X-ray photon index 
and the QPO frequencies is exp ected. Such correlation is in- 
deed observed in different XRB (|Titarchuk fc Shaposhnikovl 
l2005h . 

The position of the pivot can also be drastically changed 
and moved towards lower energies if the seed photons typi- 
cal energy is far lower than considered here . For example, as 
suggested e.g. in lMalzac fc Belmont! (|2009l ), the soft photon 
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contribution could originate from the synchrotron emission 
in optical and UV from the hot plasma. 

We also demonstrated that Cu n increases after the 
pivot whereas C qua d is decreasing again around 20 keV. 
The general increase of 8L/L with the energy after the 
pivot (visible e.g. in the rms-spectra plotted in Fig. [7} 
could be responsible in part for the observed increase of the 
variability with energy f or the QPO and the continuum (see 
e.g. for XTE J1550-564 ICui et al.|[l999h . Alternatively, the 
obse rved decrease of the QPO harmonic strength in Fig. 
5 of ICui et al.1 (|l999h could also be well explained by the 
drop in the value of C qU ad after 30 keV. We can also note 
the simila r ity be tween the RMS spectra plotted in Fig. 5 of 
ICui et al.1 (| 19991) . when the source is in its HIMS, and the 
one we obtained with our model. The increase of the QPO 
and continuum variability between 2 and 20 keV is then a 
direct consequence of the comptonisation process. 

We also observe a decrease of the powerlaw slope in the 
f-spectra when the frequencies get higher. This is consistent 
with t he results obtained for Cyg X-l in iRevnivtsev et alj 
(|l999h . It is worth noting that those authors explain the 
presence of "wiggles" in the f-spectra as a consequence 
of reflexion features on an optically thick material. Our 
Monte-Carlo simulations shows however that those features 
can be generated by the comptonisation process only. 

The non-linear domain is then investigated as a func- 
tion of the energy. We demonstrated that if a reasonable 
arbitrary value of the non-linear to linear ratio is chosen, 
it is expected that the non-linear effects of the radiative 
response will occur only at very high and medium energies. 

The choice of a cylindrical geometry has also an 
impact on the degree of loss of the input signal due to 
scatterings. The following arguments are si milar to thos e 
discussed in previous studies, such as in e.g iMiller! (|l995f ). 
Indeed, unless one global oscillation takes place in the 
corona (which is the case at very low exciting frequencies), 
a photon may encounter alternatively positive and negative 
perturbations in pressure if it travels in the radial direction. 
It would consequently kill the effect of the input signal 
on the output lightcurve. We even expect this effect to be 
emphasised when looking at high energies, as, regardless 
of the optical depth of the medium chosen, a high energy 
photon encountered a large number of scatterings. On the 
contrary, the oscillation will not be smeared out if the travel 
of the photon remains local. 

As we chose a cylindrical geometry, the optical depth 
r ~ 1 is relative to the vertical direction and hence this 
implies that in the azimuthal and radial directions r r ^ 3> 1. 
As a consequence, this is forcing the locality, relative to 
the radial wavelength of the pertubation, which in turn 
preserves the signal, whatever the output photon energy. 
As one goes towards a more spherical geometry, where r 
now becomes a measure of the radial optical depth, r in 
various directions becomes much less anisotropic. Looking 
at higher energies means looking at a larger range of 
sampled locations, and hence smearing of the imparted 
variability signal for sufficiently small wavelengths of the 
perturbation. 



explain the main features observed in the PSD of XRB, 
significant improvements are needed to include different, 
and potentially important, physical effects. 

For instance we have shown that the effect of trans- 
mitting the wave in an inner medium results in softening 
the resonances and the gaps in the power spectra. But 
the kinematics of the corona could also play an important 
role. Assuming the corona located in the inner region of an 
accretion flow, some (differential) rotation is expected and 
we might expect the emitted frequency by each part of the 
corona to be Doppler shifted, smoothing the PSD in the 
same way emission lines can be smoothed in energy spectra. 
In the close vicinity of the central engine, the general 
relativistic effects should also significantly contribute to 
blurred the PSD. The use of a non uniform sound velocity 
profile (due to e.g. non uniform corona temperature) 
along the radial or vertical directions is also expected to 
have some impact in this respect. But more importantly, 
including all these effects will result in a very different 
dispersion equation followed by the acoustic waves and 
thus to a significantly different behavior of the radiative 
response of the corona. 

We also used a single sound wave, but other acoustic 
waves could be present. For example, in the case of magnetic 
plasma Alfven and both slow and fast magnetosonic waves, 
with velocities va, c~ and c + , should be used. The former 
do not generate density or pressure perturbation. They 
are not expected to change the comptonisation efficiency, 
contrary to the magnetosonic ones. In the case of plasma in 
cquipartition, c 2 s = v\, hence c~ = 0.54c s and c + = 1.31c s 
and we do not expect significant changes compared to the 
sound wave case studied her^E However, for plasma far 
from equipartition, ~ c s ^ c + ~ va- Therefore if both 
wave propagates in the medium, each component would be 
responsible for its own band limited noise and resonances 
in the power spectra. Multiple perturbating waves may 
then explain the multiple Lorentzian components g enerally 
needed to fit the PSD o f XRB in hard states (e.g. iNowakl 
|2000| ; iPottschmidt et al.1 12003T) . Note however that in the 
HIMS, only two Lorentzians are necessary. If we identify 
in our framework t he low frequency Lorentzian (Lt,, see 
iBelloni et ail (|2002T ) for the labels of the different PSD 
components) with the propagation of a slow magnetosonic 
wave and the lower upper frequency (Li) with the fast one, 
it would require a very high magnetisation parameter for 
the plasma. Indeed the obse rvations gives rough ly vi ~ 50^6 
(see e.g. tables 2 and 3 in IBelloni et alj|2002h . and since 
c~ ~ c B and c + ~ va, it would imply that v\lc\ ~ 2500. 

We also did not investigate the implied time-lags gen- 
erated by the comptonisation process. It is however known 
that such la gs can be reproduced by piv oting in the spec- 
tra ( see e.g. iPoutanen fc Fabian|[l999l or iKording fc Falckd 
l2004h . and our models predicts the presence of such a pivot. 

More fundamentally, as we probe the timing response of 
the corona, it can be used as a complement to other variabil- 
ity models as the shape of the input excitation is supposed 
here to be a white noise. T he RMS-flux scaling foun d in 
both AGN and XRB (see e.g. lUttlev fc McHardvllioOll ) has 



While this simple model already provides very in- 
teresting timing behavior and appears very promising to 



3 The above calculations were done with va, poloidal = l/%/2 Va- 
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demonstrated that the proce sses at the origin of the variabil- 
ity is, by essence, non-linear (|Uttlev et al.ll2005h . This is not 
necessarily in contradiction with our study since we only in- 
vestigate what is the timing response of the corona, without 
taking into account the feedback process of the corona on the 
disc itself. Our model just examine the subsequent effects of 
the corona and in particular how the radiative transfer and 
the geometry of this optically thin medium could filter any 
input variability. 

The present model does not intend to reproduce all the 
observed features concerning XRB variability, especially due 
to its simpleness. The results obtained in this paper are how- 
ever very promising for further investigations, and a more 
detailed analysis, including the effect of the corona rotation 
and geometry will be presented in a forthcoming paper. 
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Using the four equivalents IA31 IA4I IA5I and IA6I the numer- 
ator N and denominator D writes: 
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Those equations can be rewritten as: 
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Finally, we get: 
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These last equations combined with Eq. 
written in Eq. [20] 
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lead to the result 



APPENDIX A: LOW FREQUENCY 
DEPENDENCE OF THE MODULATING 
FUNCTION M 

Hankel's functions can be expressed in function of Bessel's 
one: 



Hn(x) = Jo(x)-iY (x) 
H 2 n {x) = J {x) + iY (x) 



(Al) 
(A2) 



For low values of x, the equivalents of the Bessel's function 
are the following: 

Jo(x) ~ 1 (A3) 
Y (x) ~ |(M(|) +7 )J„(,) (A4) 
Ji{x) ~ | (A5) 
Yi(x) ~ -A + ^ ln (|)_^(_2 7 + l) (A6) 

TTX IT \ Z I ATT 

|Abramowitz & Stegunl Il964l ). where 7 is the Euler- 
Mascheroni constant. If the Hankel's functions are expressed 
in terms of Bessel's functions, both numerator and denom- 
inator of the function M (Eq. 1 1 8 p are pure imaginary and 
the modulating function M writes: 

„ Jo,i(rij-en,i)-iWij-e/i,i) 
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